CANONICAL FORMS OF 2x2x2 AND 2x2x2x2 SYMMETRIC 
TENSORS OVER PRIME FIELDS 



STAVROS G. STAVROU 

Abstract. We consider symmetric tensors of format 2x2x2 and 2x2x2x2 
over prime fields. Using computer algebra we compute the canonical forms of 
these tensors. For 2x2x2 symmetric tensors, we consider the prime fields F p 
for p = 2,3, 5, 7, 11, 13, 17. For 2x2x2x2 symmetric tensors, we consider the 
prime fields F p for p = 2, 3, 5, 7. For each canonical form, we determine the 
size of its orbit and the symmetric rank of the symmetric tensors in its orbit. 



1. Introduction 

A tensor is a multidimensional array of numbers. Formally, an order-fc tensor is 
an element of the tensor product of k vector spaces. A first order tensor is a vector 
and a second order tensor is a matrix. These objects were studied by Cayley in 
the late 1800's and remain an active area of research [7J [41]. The notion of tensor 
rank and the deomposition of a tensor into a sum of outer products of vectors was 
first introduced by Hitchcock [23] [24] . Models of this decomposition were formed 
later: Harshman named his model PARAFAC (short for parallel factor analysis) 
[2"2"] ; and Carrol and Chang named their model CANDECOMP (short for canonical 
decomposition) 6 . Now these models are collectively referred to as CPD (short for 
canonical polyadic decomposition) [27] . For a textbook on tensor decompositions 
and related algorithms, as well as applications of tensors in areas such as data 
mining, email surveillance, gene expression classification, and signal processing, see 
Cichocki et al. f9"]. Kroonenberg |28j. and Smilde et al. [39] , Symmetric tensors have 
applications in electrical engineering, signal processing, computational complexity, 
telecommunications, and data analysis; see [1] [10] [IT] [30] [31] and references 
therein. Bernardi et al. discuss algorithms for computing the symmetric rank of 
2 x • • • x 2 tensors from a geometric point of view [1] . Landsberg et al. consider 
the rank and border rank of symmetric tensors using geometric methods |31] . 

In this paper we use computer algebra to study symmetric tensors of format 
2x2x2 and 2x2x2x2 over prime fields, where we compute the canonical forms 
with respect to the action of the direct product of general linear groups GL2 ! 2.2(I r p) 
and GL 2j 2,2,2(Fp), respectively. 

More generally, the canonical forms of 2 x 2 x 2 tensors have been studied many 
times in the literature [37] [35|. In 1881, Le Paige was the first to consider this 
problem over C [32]. Over R, the problem was treated by Oldenburger [36]. More 
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recently, de Silva and Lim summarized these results |13j . In recent work, we con- 
sidered the problem of determining canonical forms of 2 x 2 x 2 and 2x2x2x2 
tensors over F2 and F3 [3J . 

For the 2x2x2x2 case, the problem of determining canonical forms over R and 
C has not yet been determined. This problem is related to quantum entanglement, 
where the classification of entanglements of four qubits corresponds to the determi- 
nation of canonical forms of these tensors. See Luque and Thibon [34] , Verstraete et 
al. [12], Borsten et al. [2], Briand et al. [I], Chterental and Djokovic [8 , Brylinski 
[5], Gour and Wallach [2D], Wallach [33] , Lamata et. al [2J5], and references therein. 

For 2x2x2 and 2x2x2x2 symmetric tensors over R and C, the problem was 
solved by Gurevich [2 lj . More recently, Weinberg [4? independently reported the 
results again. In this paper, we extend our results in [3J by considering the same 
problem, but restricting our attention to symmetric tensors. 

2. Preliminaries 

Definition 1. An order-A; tensor X is an element of the tensor product of k 

vector spaces Vi (g> V2 <S> • • • <8> Vk , where the order refers to the number of dimensions 
or modes. 

This is an abstract definition that does not depend on a choice of basis in each 
vector space. This is one of the motivations for studying the canonical forms of 
tensors, since they represent properties of tensors which do not depend on the choice 
of basis. Once we fix a basis in each vector space, we can associate to each order-Ac 
tensor a fc-dimensional array. Then an order-1 tensor is a vector, and an order-2 
tensor is a matrix, and we say that the order-Ac tensor of format d\ x • • • x dk is an 
element of x ••• xd fc _ j n particular, we consider order-A; tensors X — [x^...^] of 
format 2 x • • • x 2 (with k factors, for k = 3,4) with entries in the prime field F p , 
where i\, . . . ,i k E {1, 2}. 

Let's discuss different ways to represent tensors. One method is to write the 
entries of a tensor into vector (flattened) form. 

Definition 2. [27] The vectorization or flattening of X, denoted Vec(X) writes 
the columns of X in a vector format, where the entries are in lexical order of the 
Ac-tuples of the subscripts. 

We can also represent order-Ac, k > 2 tensors X in terms of matrix slices. We 
consider 2x2x2 tensors. 

Definition 3. A slice of 1" 6 F 2x2x2 is a 2 x 2 matrix obtained by fixing one of 
the three indices. There are two slices in each of the three modes: 

XU2 X122 
X212 X 2 22\ ' 

X\2l X122 
X221 X 2 22\ ' 

X211 X212 
X221 X 2 22_ 

We conventionally represent X E F 2x2x2 in terms of its frontal slices: 

33112 2:122 
^212 2:222 



frontal slices (fix A:): 



2:111 2:121 
2:211 2:221 



vertical slices (fix j): 



xm X112 
X211 X212 



horizontal slices (fix i) 



xm X112 
2:121 2:122 



X = 



Xui X121 

X211 X221 
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Although we make use of representing 2x2x2 tensors in terms of slices and its 
vectorized form, we restrict ourselves solely to the vectorized form of 2x2x2x2 
tensors, since it's the simplest way to visualize them. 



Definition 4. [12 An order-/c tensor [x^...^] G F r ' 



is called symmetric if 



h,...,i k G {l,...,n}, 



for all permutations 7r G Sk- 

Proposition 5 (Proposition 3.7 12 ). Let X = [x n ... lfc ] G F" x " 
tensor. Then 

tt(X) = X 

for all permutations ir G Sk if and only if 

■^V(l) • = x il---ik ' *1> • • • i ^fe S {1, . . . , 71 } 

/or aZ/ permutations n £ Sk- 

For example, an order-3 tensor [xjjfc] G F" xnx ™ is symmetric if 

X%jk •Eikj %jik — Xjki — Xkij — Xkji 

for all k G {1, . . . , n}. In particular, consider the case when n = 2 



be an order-k 



X 



Xm X-L21 
X211 X221 



X112 X122 
X212 X 2 22 



Then x\\2 — x\2\ — X211 —'■ b and #122 = 2:212 = X221 =: c, and so we can label X 
in the following way: 



X 



a b 
b c 



b c 
c d 



We will denote the set of all order- k , n-dimensional symmetric tensors over 
the field F by S k (¥ n ) C F" x ' " xn . The set of such tensors satisfies the property 
n(X) = X for all tt G S k and X G S k (¥ n ). 



is decomposable if it can be written 



Definition 6. 23] A tensor X e F dlX ' xdfc 
as 

X = (g>it (2) (g) • • • u (k \ 
with non-zero ifW G F di for i = 1, . . . , k. The 12, ... , ifc)th entry of X is 



(1) (2) (fe) 



A decomposable tensor is also called a simple or rank-1 tensor. If we let 
— . . . — u( k ) in the definition above then we have defined a symmetric 
simple tensor. This allows us to define two different notions of tensor rank. 

Definition 7. 23J A tensor has outer product rank r if it can be written as a 
sum of r (and no fewer) decomposable tensors, 



X = Yj = u i ® ' ' • ® «i = u i ® ■ • • ® 4 + • • • + 4 1} ® 



(8) u 



(fe) 



4 1} G F<V.., U < fc) 
outer product rank of X . 



where uY 1 G F rfl , . . . ,u\ G ¥ dk ,i — 1, ...,r. We write rank(X) to denote the 
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Definition 8. A tensor has symmetric outer product rank s if it can be written 
as a sum of s (and no fewer) symmetric simple tensors, 

X = ±uf k 

i=l 

We write rankg(X) to denote the symmetric outer product rank of X. 

The only rank-0 (symmetric or not) tensor is the zero tensor. We will drop the 
words outer product, and simply say rank and symmetric rank. 

Definition 9. The maximum rank is defined to be 

max{rankpO | X G Y dlX "' xdk }. 

If we replace rank(AT) with ranks(JT) we get the analogous definition of maximum 
symmetric rank. 

Given a symmetric tensor X we can compute its decomposition in S k (V n ) or in 
F™ x "' xn . Since the decomposition in the former set is constrained, it follows that 
for all X £ S k (¥ n ) we have the following inequality relating the two notions of 
rank: 

rank(X) < ranks (X). 

In order to compute the canonical forms of symmetric tensors we use a tensor- 
matrix multiplication. Multilinear matrix multiplication is a tensor-matrix 
multiplication that allows us to multiply matrices on each of the modes of the 
tensor. If X = [x h ... ik ] € Y dlX - xdk and 

Ai = [a\f]e¥ a ^ d \...,A k = [<4 fc) ] eF"' x4 , 

then Y = (Ai,. . . ,A k ) ■ X = [y tl ...i k ] 6F° lX '" >!o * is the new tensor defined by 

di,—,dk 

Uii...i k — y ] a iiji ■ ■ ■ a ikjk x jl---jk ■ 
»l «fc=l 

That is, the matrices Ai act by changes of basis along the k modes of the tensor. 
Since we are studying symmetric tensors, we impose the condition that A\ = ■ ■ ■ = 
Ak, otherwise, multilinear matrix multiplication can transform a symmetric tensor 
X into a non-symmetric tensor Y. 

The set of n x n invertible matrices with entries in F forms a group under matrix 
multiplication called the general linear group, denoted by GL„(F). 

Definition 10. The group GL du ... )dk (¥) := GL dl (¥) x • • ■ x GL dk (¥) acts on an 
order- k tensor by changes of basis along the k modes. 

We consider the action of the symmetry group GL 2j ... ; 2(F p ), which does not 
change the rank of a tensor [T3] . Moreover, the action of this group decomposes the 
set of ovdev-k symmetric tensors into a disjoint union of orbits, where the tensors 
in each orbit are equivalent under the group action. The orbit of X is the set 

Ox :={(s,s,s)-X| 5 eGL 2 (F p )}. 

We define the canonical form of X to be the minimal element in its orbit with 
respect to the lexical ordering. 
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symmetric rank canonical form 





n n 

u u 



n n 





1 


' 1 



" 





2 


' 1 




" 
1 




3 

3 [ ° 
[ 1 


' 1 

1 

1 

-1 


1 " 



1 - 

-1 


1 " 



Table 1. Canonical forms of 2 x 2 x 2 symmetric tensors over R 



Every order-fc symmetric tensor of dimension n may be uniquely associated with 
a homogeneous polynomial of degree k in n variables [12J . Homogeneous polyno- 
mials are also called quantics. Then the problem we consider is equivalent to the 
problem of determining the minimal number of pth powers of linear terms: the Big 
Waring Problem. This has been considered by many authors, including Ehrenborg 

nana- 

3. Symmetric Tensors of Format 2x2x2 

3.1. R and C. Gurevich and Weinberg [33J computed the canonical forms of 
symmetric 2x2x2 tensors over R and C by considering them as homogeneous 
polynomials. Over R the canonical forms of X E <S 3 (R 2 ) and the symmetric rank 
are listed in Table Q] 

Over C, the list is the same except that we omit the tensor 



" 1 


1 -1 " 


1 -1 


-1 



so that there is only one orbit in symmetric rank 3. The maximum symmetric rank 
is 3 over R and C. We will see shortly that the maximum rank is larger over some 
of the prime fields ¥ p . 

3.2. Algorithms. The set of 2 x 2 x 2 tensors over the prime field F p contains p 8 
tensors: p many elements in each of the eight entries of the tensor. If we restrict 
ourselves to the symmetric case, there are only p A tensors. The algorithm has two 
main components: generate all the tensors in each rank; apply the group action so 
that the set of tensors in each rank is decomposed into a disjoint union of orbits. 
We've used this algorithm in previous work to determine the canonical forms of 
2x2x2 and 2x2x2x2 tensors over finite fields 3 j . There are two modification: 
in the first part, the set of rank-1 tensors must contain only symmetric simple 
tensors; in the second part, the same group element must act on every mode of 
the tensor whose orbit we wish to compute (otherwise the orbit will include non- 
symmetric tensors). 
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f latten(.x) 

return([a; ul) x n2 , x 121 , x 122 , x 2n , x 212 , x 221 , x 222 }) 

outerproduct(a, b, c) 

for i = 1, 2 do for j = 1, 2 do for k = 1, 2 do: set lyj. -s— CLia^a^ mod p 
return( £ ) 

• set vectors «- { [«, j] | < i, j < p-1} \ {[0,0]} 

• set arrayset[0] <- {[0,0,0,0,0,0,0,0]} 

• set arrayset[l] <— { } 

• for a in vectors do 

— set x f latten(outerproduct(a, a, a)) 

— if x arrayset[0] and x ^ arrayset[l] then 

set arrayset[l] <— arrayset[l] U {x} 

• set r <— 1 

• while arrayset[r] ^ { } do: 

— set arrayset[r+l] 4— { } 

— for x G arrayset[r] do for y 6 arrayset[l] do 

* set z <!— [ a;i+yi modp, . . . , xs+2/8 modp] 

* if z ^ arrayset[s] for s = 0, . . . , r+1 then 

set arrayset[r+l] <— arrayset[?'+l] U {z} 

— set r ^— r + 1 

• set maximumrank ^— r — 1 

Table 2. Algorithm to generate tensors in each rank (pseudocode) 



Algorithm 11. [3] 

• The only rank-0 tensor is the zero tensor. The set of symmetric rank-1 tensors is 
{a ® a g) o | a g Fy} where a is a non-zero vector. 

• Assume that the symmetric rankg-r tensors have already been computed. To 
compute the symmetric tensors of rankg-(r + 1), we consider all sums X + Y where 
rankg (A) = r and ranks (F) = 1. Clearly, ranks (A + Y) < r + 1, but it is 
possible that ranks (A + Y) < r + so we only retain those A + Y such that 
ranks(A + F) = r+l 

• Sort the tensors within each symmetric rank in lexical order. 

• For each symmetric rank, perform the following iteration until there are no more 
tensors in the current symmetric rank: 

- Choose the minimal element of the set of tensors of the current symmetric 
rank. 

- Compute the orbit of this element under the action of the symmetry group. 

- Remove the elements of this orbit from the set. 

See Table [2] and Table [3] for the pseudocode from [3] that has been modified for 
symmetric tensors. 

3.3. F2. Over F2, there are 256 tensors, 16 of which are symmetric. There are three 
symmetric simple tensors, 



" 1 








" 




" 








" 




' 1 1 


1 1 " 

























1 




1 1 


1 1 
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unf latten(x) 

set t 4- 

for i = 1, 2 do for j = 1, 2 do for fc = 1,2 do: set i •<— t + 1; set x t 
return(y) 

groupaction(5, x, to) 
set j/ ^— unf latten(x) 

if m = 1 then for j = 1, 2 do for fc = 1, 2 do: 

set v <- [yijk, V2jk}; set w <- [gnv 1 +g 12 v 2 modp, 521^1+522^2 modp] 

for i = 1, 2 do: set j/^ «— Wi 
if m = 2 then . . . (similar for second subscript) 
if m = 3 then . . . (similar for third subscript) 
return( f latten( y ) ) 

smallorbit(x) 
set result <— { } 
for a e GL 2 (¥ p ) do: 

set y <— groupaction(a, x, 1) 

set z ^— groupaction(a, y, 2) 

set w <~ groupaction(a, z, 3) 
set result result U {w} 
return( result ) 

• for r = 0, . . . , maximumrank do: 

set representatives^] -«—{}; set remaining arrayset[r] 
while remaining 7^ { } do: 

set x remaining[l]; set xorbit <— largeorbit(x) 

append xorbit[l] to representatives^] 

set remaining <— remaining \ xorbit 

Table 3. Algorithm for group action (pseudocode) 



We emphasize that not all 16 symmetric tensors can be written as a sum of sym- 
metric simple tensors over F 2 . The following 8 tensors do not have a symmetric 
decomposition: 



" 


1 


1 


" 




' 1 


1 


1 


" 




' 


1 


1 


" 




' 1 


1 


1 


" 


1 











•) 


1 











7 


1 








1 


) 


1 








1 


" 








1 " 




' 1 








1 " 




' 








1 " 




' 1 








1 " 





1 


1 





•) 





1 


1 





7 





1 


1 


1 


) 





1 


1 


1 



The maximum symmetric rank over F2 is 3. The number of symmetric tensors 
in each symmetric rank and the approximate percentages are listed below. The 
percentages add to 50%, since only half of the tensors have a decomposition into a 
sum of symmetric simple tensors. 

rank 12 3 

number 13 3 1 

w % 6.25% 18.75% 18.75% 6.25% 
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symmetric rank orbit size canonical form 



1 


' 



" 





1 3 


' 



" 
1 




2 3 


" 1 
1 1 


1 1 " 
1 1 




3 1 


' 1 

1 1 


1 1 " 

1 





Table 4. Canonical forms of 2 x 2 x 2 symmetric tensors over F2 



symmetric rank orbit size canonical form 



1 


' 



" 





1 8 


' 




" 
1 




2 24 


' 1 
1 


1 " 
1 




3 8 
3 24 


' 

1 
: 1 

1 


1 " 

1 

1 = 
2 




4 16 


' 
1 


1 " 

1 1 





Table 5. Canonical forms of 2 x 2 x 2 symmetric tensors over F3 



The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table 21 

3.4. F 3 . There are 6,561 tensors over F 3 , 81 of which are symmetric. Unlike in 
the previous case, all 81 symmetric tensors can be generated as a sum of simple 
symmetric tensors. The maximum symmetric rank in this case is 4. The number 
of symmetric tensors in each symmetric rank and the approximate percentages are 
listed below. 

rank 12 3 4 

number 1 8 24 32 16 

% 1.23% 9.88% 29.63% 39.51% 19.75% 

The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table [5j 
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symmetric rank orbit size canonical form (flattened) 






1 


[o 




















0" 


1 


24 


[o 




















1] 


2 


240 


[o 


1 


1 





1 








1] 


3 


120 











1 





1 


1 


0" 




3 


80 





1 


1 





1 








2^ 




3 


160 


1 








1 





1 


1 


2" 





Table 6. Canonical forms of 2 x 2 x 2 symmetric tensors over F5 



3.5. F5. There are 390,625 tensors over F3, 625 of which are symmetric. The max- 
imum symmetric rank in this case is 3. The number of tensors in each symmetric 
rank and the approximate percentages are listed below. 

rank 12 3 

number 1 24 240 360 

« % 0.16% 3.84% 38.40% 57.60% 

The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table [B] From now on, the canonical forms will be displayed 
in flattened form. 

3.6. F7. There are 5,764,801 tensors over F7, 2,401 of which are symmetric. The 
maximum symmetric rank in this case is 5, and we also see that the number of 
orbits has jumped to 14. The number of symmetric tensors in each symmetric rank 
and the approximate percentages are listed below. 

rank 1 2 3 4 5 

number 1 16 128 688 1232 336 

w % 0.04% 0.67% 5.33% 28.65% 51.31% 13.99% 

The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table 

3.7. Fn. There are 214,358,881 tensors over Fn, 14,641 of which are symmetric. 
The maximum symmetric rank in this case is 3. The number of symmetric tensors 
in each symmetric rank and the approximate percentages are listed below. 

rank 12 3 

number 1 120 6600 7920 

« % 0.01% 0.82% 45.08% 54.09% 

The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table [8j 

3.8. F13. There are 815,730,721 tensors over F13, 28,561 of which are symmetric. 
The maximum symmetric rank in this case is 4, and there are 14 orbits. The number 
of symmetric tensors in each symmetric rank and the approximate percentages are 
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symmetric rank 


orbit size 


canonical form (flattened) 
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1 


n 

u 


n 

u 


n 

u 


n 


n 

u 


n 
u 


n n 
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16 
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112 
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3 


16 




















3 


3 


112 





1 


1 





1 





1 


3 


336 





1 


1 





1 





3 


3 


224 


1 

















2 


4 


336 











1 





1 


1 


4 


112 





1 


1 





1 





6 


4 


336 





1 


1 





1 





6 


4 


224 


1 

















3 


4 


224 


1 








1 





1 


1 2 


5 


336 





1 


1 





1 





5 



Table 7. Canonical forms of 2 x 2 x 2 symmetric tensors over F7 



symmetric rank orbit size canonical form (flattened) 
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1 
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2 


6600 
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1 
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3 


1320 











1 





1 


1 





3 


2200 





1 


1 





1 








2 


3 


4400 


1 








1 





1 


1 


2 



Table 8. Canonical forms of 2 x 2 x 2 symmetric tensors over Fn 



listed below. 

rank 12 3 4 

number 1 56 1568 16016 10920 

w % 0.00% 0.20% 5.49% 56.08% 38.23% 

The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table GO 

3.9. F17. There are 6,975,757,441 tensors over F13, 83,521 of which are symmetric. 
The maximum symmetric rank in this case is 3, and there are 6 orbits. The number 
of symmetric tensors in each symmetric rank and the approximate percentages are 



CANONICAL FORMS OF 2x2x2 AND 2x2x2x2 SYMMETRIC TENSORS 11 



symmetric rank orbit size canonical form (flattened) 
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Table 9. Canonical forms of 2 x 2 x 2 symmetric tensors over F13 



symmetric rank orbit size canonical form (flattened) 
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Table 10. Canonical forms of 2 x 2 x 2 symmetric tensors over F17 



listed below. 

rank 12 3 

number 1 288 39168 44064 

« % 0.00% 0.34% 46.90% 52.76% 
The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table [TD] 

For p = 19 there are almost 17 billion tensors, where over 130,000 of these 
are symmetric. For this and larger primes, our program does not have sufficient 
memory. 

4. Symmetric Tensors of Format 2x2x2x2 

4.1. R and C. Weinberg [44] lists the canonical forms of 2x2x2x2 symmetric 
tensors based on the discussion from Gurevich [21] . Recall we represent these fourth 
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order tensors as X = [xijki] where i,j,k,l = 1,2. Over R there are 9 canonical 
forms: 

1. xmi = 1, 2^(1122) = M, ^2222 = 1, \i < -|, itjfci = otherwise; 

2. Xmi = a, ^0.(1122) = "M; ^2222 = «, |U > _ |) a = ±1- ^ij'fci = 

otherwise; 

3. inn = 1, a; (T (ii22) = Mj ^2222 = -1, Xijki = otherwise; 

4. x a (xxi2) = a, Xijki = otherwise; 

5- Xijki = 0; 

6- x a{1122 ) = ±1, X2222 = ±1, ^ijfc; = otherwise; 

7- 2V(ii22) = a, x^ki = otherwise; 

8. xnn = 7, a; (T ( 11 22) = £2222 = 7,7= 1, Xyjw = otherwise; 

9. Xmi = ±1, Xijki = otherwise. 

Here, a S S4 denotes any permutation of the indicated subscripts. If we consider 
the group action 01^2,2,2,2 (C) then we get 6 canonical forms: 

1. xmi = 1, ££,(1122) = 6/i, .T2222 = 1, M 7^ ^ti = otherwise; 

2- £ CT (U22) = 6, X2222 = 1, x ijk i = otherwise; 

3- £ CT (U22) = 6, x ijk i = otherwise; 

4- x er(ul2 ) = 4, iy t! = otherwise; 

5. xmi = 1, x^ki = otherwise; 

6. Xijki = 0. 

The algorithms we used for the 2x2x2 case generalize to the 2x2x2x2 case 
by making a few simple modifications. In particular, we consider the outer product 
of four non-zero vectors (instead of 3), and we add an extra index £ 6 {1,2} for the 
group action on the fourth mode. 

4.2. F2. There are 65,536 tensors over F2, 32 of which are symmetric. Just as we 
experienced for the 2x2x2 case over F2, not all the symmetric tensors can be 
represented as a sum of symmetric simple tensors. In this case, only 8 out of 32 
symmetric tensors can be represented as a sum of symmetric simple tensors. The 
maximum symmetric rank in this case is 3, and there are 4 orbits. The number 
of symmetric tensors in each symmetric rank and the approximate percentages are 
listed below. 

rank 12 3 

number 13 3 1 

w % 3.125% 9.375% 9.375% 3.125% 

The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table [TT] To improve legibility, we use a period (.) in place 
of zero (0). 

4.3. F3. There are 43,046,721 tensors over F3, 243 of which are symmetric. Only 
81 out of the 243 symmetric tensors can be written as a sum of symmetric simple 
tensors. The maximum symmetric rank in this case is 8, and there are 15 orbits. 
The number of symmetric tensors in each symmetric rank and the approximate 
percentages are listed below. 

rank 012345678 
number 1 4 10 16 19 16 10 4 1 

ra % 0.41% 1.65% 4.12% 6.58% 7.82% 6.58% 4.12% 1.65% 0.41% 
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Table 12. Canonical forms of 2x2x2x2 symmetric tensors over F3 



The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table [12] 

4.4. F5. There are 152,587,890,625 tensors over F 5 , 3,125 of which are symmetric. 
In this case, all symmetric tensors have a symmetric decomposition. The maximum 
symmetric rank in this case is 10, and there are 52 orbits. The number of symmetric 
tensors in each symmetric rank and the approximate percentages are listed below. 
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Table 13. Canonical forms of 2x2x2x2 symmetric tensors over F5 



The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Tables [13] and [TU 



4.5. F7. There are 33,232,930,569,601 tensors over F7, 16,807 of which are symmet- 
ric. Again, all symmetric tensors have a symmetric decomposition. The maximum 
symmetric rank in this case is 6, and there are 42 orbits. The number of symmetric 
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Table 14. Canonical forms of 2x2x2x2 symmetric tensors over 
F 5 (cont.) 



tensors in each symmetric rank and the approximate percentages are listed below. 

rank 1 2 3 4 5 6 

number 1 24 276 1932 7119 6615 840 

w % 0.01% 0.14% 1.64% 11.50% 42.38% 39.36% 5.00% 

The symmetric ranks, orders of each orbit, and the minimal representatives of 
each orbit are given in Table [15] 

We mention now that for larger primes our computer did not have sufficient 
memory to complete the computations in a reasonable amount of time. 

5. Conclusion 

In this paper we began by summarizing known results about the canonical forms 
of 2 x 2 x 2 and 2x2x2x2 symmetric tensors over R and C. Our original 
contribution was considering the same problem over specific prime fields ¥ p . We 
obtained a classification of canonical forms of symmetric tensors, and verified that 
the maximum symmetric rank over certain prime fields is greater than that over 
R and C. Based on the prime fields we considered, the maximum symmetric rank 
of 2 x 2 x 2 symmetric tensors is at least 5. The maximum symmetric rank of 
2x2x2x2 symmetric tensors is at least 10. Due to memory limitations, we could 
not consider larger prime fields. 

Our algorithms are modifications of those used in [3] , where we computed canon- 
ical forms of 2 x 2 x 2 and 2x2x2x2 tensors over F2 and F3. 
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Table 15. Canonical forms of 2x2x2x2 symmetric tensors over F7 
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